Self-affine surface morphology of plastically deformed metals 
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We analyze the surface morphology of metals after plastic deformation over a range of scales 
from 10 nm to 2 mm, using a combination of atomic force microscopy and scanning white-light 
interferometry. We demonstrate that an initially smooth surface during deformation develops self- 
affine roughness over almost four orders of magnitude in scale. The Hurst exponent H of one- 
dimensional surface profiles is initially found to decrease with increasing strain and then stabilizes 
at H » 0.75. By analyzing their statistical properties we show that the one-dimensional surface 
profiles can be mathematically modelled as graphs of a fractional Brownian motion. Our findings 
can be understood in terms of a fractal distribution of plastic strain within the deformed samples. 



According to the traditional paradigm implicit in con- 
tinuum models of plasticity, plastic deformation under 
homogeneous loads and in the absence of so-called plas- 
tic instabilities is expected to proceed in a smooth and 
spatially homogeneous manner. Fluctuations, if any, are 
supposed to average out above the scale of a "representa- 
tive volume element" which is tacitly assumed to be small 
in comparison with macroscopic dimensions (or, indeed, 
in comparison with the scale of any problem to which 
continuum modelling is applied). Recently for crystalline 
systems deforming by dislocation glide this paradigm has 
been challenged both from the theoretical and experi- 
mental side 0, 0> 0, 0> IS El • Weiss et. al. demonstrated 
that dislocation glide in single-glide oriented ice single 
crystals proceeds in a temporally intermittent manner, 
i.e., plastic flow is characterized by a sequence of large 
"deformation bursts" 0, 0]. Acoustic emission due to 
these bursts was monitored and it was demonstrated that 
there is no characteristic burst magnitude. Instead the 
energy releases were found to obey a scale- free power-law 
distribution over more than six orders of magnitude in 
energy 0] . Simulations of dislocation dynamics in single- 
glide single crystals showed qualitatively similar behavior 
mm These findings were interpreted as dynamical criti- 
calphenomena in a slowly driven non-equilibrium system 
0, , and it has been argued that the dislocation system 
in a plastically deforming crystal is always close to a crit- 
ical point ("yielding transition", (4(). 

Temporal intermittency of plastic flow goes along with 
spatial heterogeneity. This has been long known to met- 
allurgists who observed pronounced surface traces ("slip 
lines" ) resulting from the collective motion of many dislo- 
cations. By monitoring the evolution of surface features 
it is therefore possible to study collective behavior in plas- 
tic flow. However, most investigations (for an overview, 
see 7J) have focused on the evolution of single slip lines 
or slip line bundles, i.e., the motion of isolated dislocation 
groups. Long-range correlations within the slip pattern 
have only occasionally been studied. Konstantinidis and 
Aifantis [8j used a wavelet transform to investigate the 
slip pattern of a Fe3 Al single crystal on multiple scales be- 
tween 0.1 and 50 micrometers but did not quantitatively 



analyze their data in view of scale-invariant behavior. 
Kleiser and Bocek [j| report a self-similar pattern of slip 
lines in Cu single crystals on scales between 0.06 and 2 
micrometers and determine a fractal dimension Dp ~ 0.5 
for the set of intersection points between the slip lines and 
a line normal to the slip direction. 

Recently, Weiss and Marsan |(J investigated the three- 
dimensional patterning of slip on macroscopic scales. By 
recording acoustic emission with multiple transducers 
and using spatial triangulation, they determined the 3D 
distribution of slip bursts in the volume of an ice single 
crystal. This investigation demonstrated that the slip 
burst pattern on macroscopic scales (0.6 mm to 10 mm) 
exhibits features of a fractal set with correlation dimen- 
sion Dq ~ 2.5. However, due to the limited spatial res- 
olution of the technique no information about the slip 
distribution on smaller scales could be obtained. 

In the present study we pursue a slightly different ap- 
proach for assessing spatial correlations in plastic defor- 
mation from the nanoscopic to the macroscopic scale. 
Instead of directly studying the distribution of plastic 
strain in the bulk or on the surface of the sample, we 
monitor the evolution of the surface profile. Provided 
that the initial surface is sufficiently smooth, the evolu- 
tion of a one-dimensional profile can be directly related 
to the plastic distortion at the surface: If we define a 
local coordinate system such that the x direction corre- 
sponds to the direction of the profile and the y direction 
to the surface normal, then the derivative y x = dy/dx of 
the profile y{x) equals the component j3 yx of the plastic 
distortion tensor. Long-range correlations in the plastic 
strain pattern can be detected since they give rise to char- 
acteristic modifications of the surface morphology - more 
specifically speaking, power-law correlations in the spa- 
tial distribution of plastic strain lead to the emergence of 
self-affine surface profiles with a Hurst exponent H > 0.5 
as will be discussed below. 

A self-affine profile y{x) is characterized by the prop- 
erty that the statistical properties of the profile re- 
main unchanged under the scaling transformation x — > 
Ax, y — > X H y. The roughness exponent or Hurst expo- 
nent H can be related to a fractal dimension Dp ("box 
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FIG. 1: Stress-strain curve of a poly crystalline Cu specimen; 
deformation at room temperature, imposed deformation rate 
during loading e = 1.5 x 10 -3 s _1 ; surface profiles were taken 
before deformation as well as after unloading at the points 
indicated by the open circles. 



dimension") of the profile through Dp = 2 — H [Io|. 
Self-affine behavior is abundant in natural surfaces, a 
most prominent example in materials physics are frac- 
ture surfaces whose large-scale roughness can be charac- 
terized by a universal (material-independent) exponent 
H w 0.8 [ll|. In the present letter we report for the 
first time an analysis of the surface structure of plasti- 
cally deformed metal samples in terms of self-affine prop- 
erties. We applied a combination of atomic force mi- 
croscopy (AFM) and scanning white-light interferometry 
(SWLI) to quantitatively characterize the surface mor- 
phology over a range of scales between 10 nm and 2 mm. 
Polycrystalline copper samples of 99.9 % nominal purity 
and an average grain size of 40 /an were electropolishcd to 
a typical rms roughness of 44 nm as determined by SWLI 
over a surface area of 40/mi x40/im. Two complementary 
methods were used for surface imaging: (i) AFM scans 
(512 x 512 pixels) were taken over areas of 26^m x26^m 
and 6/^m x6/^m using a PicoSPM (Molecular Imaging) 
at the constant deflection mode with a medium range 
scanner and a cantilever (Mikromasch, Ultrasharp) with 
a nominal spring constant of 1.75 N/m and tip curvature 
radius 10 nm; the rms vertical noise amplitude was less 
than 0.05 nm; (ii) SWLI (New View 100, Zygo) scans 
were performed over areas of 130/im x2000^m with a 
lateral resolution of 0.6 fim and a vertical resolution of 
0.5 nm. From the scans, typically 4-5 surface profiles 
were taken in the direction of the specimen axis and 2-3 
profiles normal to that direction. 

After characterizing the initial surface, the samples 
were deformed in tension on a standard tensile-testing 
machine (Instron series 3360). Testing was done at 
room temperature with an imposed strain rate e = 
1.5 x 10~ 3 s _1 . At total (engineering) strains e — 
2.3%, 5.5%, 9.6%, 17.8% and 23% (see Figured, the sam- 
ple was unloaded and AFM and SWLI scans were taken 
in a similar manner as before deformation. The sequence 
of surface morphology investigations was ended at the 
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FIG. 2: Surface profiles taken at e = 9.6%; top: AFM pro- 
file; bottom: SWLI profile; the x direction is parallel to the 
direction of the tensile axis; roughness plots corresponding to 
these profiles are shown in Figure 

onset of macroscopic deformation localization (necking) . 

Typical surface profiles are shown in Figure [3 rep- 
resenting a "long" AFM profile obtained from a 26/im 
x26/im scan together with a SWLI profile. To investi- 
gate a possible self-affine structure, mean height differ- 
ences (\y(x) — y(x+L)\) were evaluated as a function of L 
for each of the AFM and SWLI profiles by averaging over 
all pairs [x, x + L) within the profile. Self-affine behavior 
implies that 

(\y(x)-y(x + L)\) K I fl . (1) 

Hence, double-logarithmic plots of (\y(x) — y(x + L)\) 
vs. L should exhibit a linear scaling regime with slope 
H. Such 'roughness plots' are shown in Figure Typ- 
ical plots exhibit linear scaling regimes extending be- 
tween 0.05 and 5 /im for the AFM and between 0.5 
and 100 (im for the SWLI profiles. Hurst exponents de- 
duced from the slope of the scaling regimes are similar 
for AFM and SWLI profiles, and also the absolute values 
of (\y(x) — y(x + L)\) are similar across the region of over- 
lap of the roughness plots. Hence, we observe continuous 
scaling over almost four orders of magnitude. Surface 
profiles y(z) taken in the direction normal to the ten- 



3 



sile axis exhibit scaling with similar exponents, although 
the length of the self-affine scaling regimes is reduced. 
Roughness plots obtained from the "short" AFM profiles 
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FIG. 3: Roughness plots (mean height difference vs. distance 
along profile) for AFM and SWLI profiles obtained at strains 
of 9.6 and 17.8 %; the corresponding profiles for e = 9.6% are 
shown in Figure 

deduced from 6/im x6/zm scans show a crossover to a 
strain-independent apparent Hurst exponent ff w 1 be- 
low 50-100 nm. 

Hurst exponents H and fractal dimensions Dp = 2 — 
H determined from the roughness plots are compiled in 
Figure 0| The error bars reflect the scatter of exponents 
obtained from different profiles taken at the same strain. 
It can be seen that the exponents initially decrease with 
increasing strain and then stabilize above e w 10% at a 
value of H f» 0.75. It is important to note that the initial 
value of H rs 0.97 does not reflect any self-affine surface 
roughness but simply stems from large-scale gradients 
on an almost smooth surface. The same is true for the 
steepening of roughness plots obtained from the "short" 
AFM profiles on very small scales: Below 50 nm, the 
apparent Hurst exponent increases towards H ~ 1 as one 
essentially observes the smooth initial surface in between 
the deformation-induced surface steps. 

To further investigate the statistical surface proper- 
ties in the different scale regimes, probability distribu- 
tions p(AyL) of surface height differences Ay^ — y(x) — 
y(x + L) have been determined for the profiles taken at 
e = 9.6%. Lq — 0.5/^m was taken as a reference and 
all height differences Ay^ were normalized by the vari- 
ance yo = (Ay£ ) ' 5 = 27.08 nm of the corresponding 
height difference distribution. Within the scaling regime 
£o < L < £i (£o w 50 nm, £i rj 100/im), probability dis- 
tributions corresponding to different values of L can be 
collapsed by re-scaling AyL/yo — * (L/Lq) h AyL/yo with 
H = 0.72 (see FigureEJ: 



p([L/L ] H Ay L /y ) =p(Ay L Jy Q ) 



(2) 
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FIG. 4: Fractal dimension and Hurst exponent as a function 
of strain, squares: values determined from scanning interfer- 
ometry profiles, circles: values determined from AFM profiles; 
each value represents an average over 5 profiles; full line: guide 
to the eye. 
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To assess the shape of the probability distributions, cu- 
mulative distributions normalized to unit variance have 
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FIG. 5: Collapse of probability distributions of height dif- 
ferences after re-scaling; data points: distributions for dif- 
ferent values of L obtained from AFM and SWLI profiles 
at e = 9.6%; insert: distributions (normalized to unit vari- 
ance) for L = 40nm (dashed line), L = 5fim (full line) and 
L — 160/^m (dotted line). 



been represented on a probability scale (i.e., the error 
function appears as a straight line). The insert of Figure 
[S] shows the cumulative distribution for L = 2/xm (full 
line) as well as two distributions for L values outwith 
the scaling regime, namely for L — 40nm (dashed line) 
and L = 160/im (dotted line). Whereas distributions 
throughout the self-affine scaling regime are well approx- 
imated by Gaussian statistics, the tails of distributions 
outwith the scaling regime exhibit characteristic devia- 
tions from Gaussian behavior. On scales below £o the 
surface morphology is characterized by flat regions sep- 
arated by huge slip steps at locations where slip events 
have occurred. This leads to a leptokurtic distribution 
with heavy tails. On scales above £i, on the other hand, 
the surface flattens and large height fluctuations are sup- 
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pressed since the macroscopic slope of the SWLI profiles 
is by definition set to zero. This gives rise to platykurtic 
distributions with truncated tails. 

The fact that the statistics is Gaussian throughout 
the scaling regime implies that the surface profiles y{x) 
can be mathematically interpreted as graphs of fractional 
Brownian motions. This allows us to draw conclusions on 
the underlying plastic strain pattern since the increments 
y x of a fractional Brownian motion exhibit long-range 
correlations given by ^3] 

{y x {x)y x (x'))^\x-x'\ 2H - 2 . (3) 

Since y x equals the component /3 yx of the plastic distor- 
tion tensor at the surface, our observations yield strong 
evidence for power-law correlations in the distribution of 
plastic strain. From the correlation integral 

I(R) = I {y x {x)y x {x'))A{x - x 1 ) oc R 2 "- 1 (4) 

J\x—x'\<R 

we deduce a correlation dimension Dq = 2H—1 ~ 0.5 for 
the one-dimensional strain pattern underlying the sur- 
face profile. The value Dq ~ 0.5 is consistent with 
the fractal dimension of slip line patterns reported by 
Kleiser and Bocek A dimension of 0.5 for a one- 
dimensional strain pattern is also consistent with the di- 
mension Dp — 2.5 for the three-dimensional pattern of 
slip events measured by Weiss and Marsan These 
investigations were performed on single crystals; due to 
experimental limitations the scaling regimes cover only 
about one order of magnitude. Our results demonstrate 
that the strain pattern in polycrystals exhibits an anal- 
ogous self-similar behavior which by combining different 
experimental techniques can be traced over almost four 
orders of magnitude in scale. 

More investigations are required to decide whether the 
Hurst exponent H « 0.75 which we determine for de- 
formed metal surfaces is universal, i.e. independent of 
material and deformation conditions. Further investi- 
gations are also needed to establish the nature of the 
boundaries of the self-affine scaling regime. The lower 
bound of the scaling regime can be related to the spac- 
ing between individual 'slip events' between which the 



surface remains unchanged by the deformation. The up- 
per boundary £i of the scaling regime, which defines the 
correlation length of the self-affine surface pattern, is of 
the same order of magnitude as the grain size (40 mi- 
crons). One may reasonably conjecture that the range of 
correlations in the dislocation dynamics, and therefore in 
the slip pattern, is limited by the grain size. However, 
further experimental work (preferably on single crystals) 
is required to decide whether this is indeed the physical 
mechanism which determines the correlation length. 

From the theoretical side, the fact that the Hurst expo- 
nent of deformed metal surfaces is quite similar to that of 
fracture surfaces raises the fundamental question whether 
there is a universal mechanism underlying both observa- 
tions, or whether the similarity of exponents is incidental. 
Crystal plasticity is a bulk phenomenon which only indi- 
rectly affects the surface, whereas fracture is directly gov- 
erned by surface and interface properties. However, this 
does not preclude the possibility of a common theoret- 
ical treatment within the framework of non-equilibrium 
critical phenomena. 

In view of plasticity models using coarse-grained strain 
variables, our results indicate that strain fluctuations 
may average out rather slowly. The strain fluctuation on 
scale L can be estimated as Aez, ~ Ay^/L. At the lower 
end of the self-affine scaling regime (£q ~ 50 nm) , fluctu- 
ations are of the same order of magnitude as the average 
strain. Within the scaling regime, strain fluctuations de- 
crease like ((,0/L) 1 ~ H if averaged over the length L; for 
(1 — H) 0.25 we find that the residual strain fluctua- 
tions at the upper end of the scaling regime (£1 « 100/im) 
are still about 15 % of the average strain. This observa- 
tion explains why it is possible to detect significant strain 
fluctuations even on macroscopic scales [Hill. 
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